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Introducing nanoparticles into superconducting materials has emerged as an efficient route to 
enhance their current-carrying capability. We address the problem of optimizing vortex pinning 
landscape for randomly distributed metallic spherical inclusions using large-scale numerical simu¬ 
lations of time-dependent Ginzburg-Landau equations. We found the size and density of particles 
for which the highest critical current is realized in a fixed magnetic field. For each particle size and 
magnetic field, the critical current reaches a maximum value at a certain particle density, which 
typically corresponds to 15-23% of the total volume being replaced by nonsuperconducting mate¬ 
rial. For fixed diameter, this optimal particle density increases with the magnetic held. Moreover, 
we found that the optimal particle diameter slowly decreases with the magnetic held from 4.5 to 
2.5 coherence lengths at a given temperature. This result shows that pinning landscapes have to be 
designed for specihc applications taking into account relevant magnetic held scales. 

PACS numbers: 74.20.De, 74.25.Sv, 74.25.Wx 


High-current applications of superconductivity 
strongly depend on engineering defect microstructures 
which efficiently suppress the mobility of vortex lines 
over a wide range of magnetic fields. Introducing 
self-assembled inclusions into cuprate superconducting 
materials has been established as a very efficient route to 
improve their current-carrying capacity. In particular, 
such inclusions may be prepared in the form of almost 
spherical particles [1-6], nanorods [7], or combinations 
of both [8]. This technology has been implemented in 
the second-generation superconducting cables based on 
YBa 2 Cu 307 (YBCO) coated conductors [9]. 

At present, improvement of pinning properties is 
achieved mostly by trial-and-error approach. Ratio¬ 
nal design of optimal pinning landscapes may be facil¬ 
itated by large-scale numerical simulations. Real strong- 
pinning materials always contain many different defects, 
such as dislocations, twin boundaries, stacking faults, etc. 
It is very challenging to understand details of vortex dy¬ 
namics in such complicated pinning landscape. While the 
ultimate critical-current optimization can be achieved by 
the constructive combination of different pinning centers, 
a natural first step is to determine the best pinning con¬ 
figuration for a relatively simple system with only one 
type of defects. In this article, we explore the case of 
monodisperse spherical defects with the size of a few 
coherence lengths which model nanoparticles in coated 
YBCO conductors. In addition, similar pinning centers 
in the form of impurity clusters can be introduced by pro¬ 
ton irradiation which provides an opportunity for further 
enhancement of the critical current [10, 11]. 

The quantitative description of vortex dynamics is a 
highly nontrivial problem even for a model system. Pin¬ 
ning is a collective phenomenon controlled by the interac¬ 


tion of vortices with pinning centers as well as the flexibil¬ 
ity of vortex lines and intervortex interactions. An ana¬ 
lytical treatment of this problem is limited to qualitative 
estimates of critical currents in simple situations, such as 
weak pinning by high densities of atomic impurities [12] 
or strong pinning by low densities of strong inclusions 
interacting with a regular vortex lattice [13, 14], see also 
reviews [15-18]. While providing a useful general guid¬ 
ance, qualitative estimates inevitably rely on simplifying 
assumptions and miss important details. 

To improve a quantitative description of pinning mech¬ 
anisms, vortex dynamics is extensively explored by nu¬ 
merical simulations using several models with various de¬ 
grees of complexity and realisticity. The minimal ap¬ 
proach is to treat the vortices as elastic strings whose 
dynamics is described by the overdamped equation of mo¬ 
tion taking into account interaction with pinning centers 
and thermal Langevin forces. Such a Langevin-dynamics 
approach [19-24] provides a qualitative description of the 
vortex state for low fields, when the distance between vor¬ 
tices is much larger than the coherence length, and for a 
low density of pinning centers. This description, however, 
has several limitations: vortex-vortex and vortex-defect 
interactions can be treated only approximately, and vor¬ 
tex cuttings and reconnections are completely neglected. 
Also, these simplified models cannot be used when pin¬ 
ning centers occupy a noticeable fraction of superconduc¬ 
tor’s volume corresponding to strongest pinning. There¬ 
fore, it is desirable to probe the strong-pinning regime by 
more sophisticated models. 

These limitations are overcome in the time-dependent 
Ginzburg-Landau (TDGL) model [25] which describes 
dynamics of the superconducting order parameter. In 
this description, the vortex lines appear spontaneously 
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Figure 1. (a) Current-voltage dependences computed for different numbers of particles A^p with diameter a = 4^ for magnetic 

field 0.1-f/c2- The numbers in parentheses represent the volume fractions occupied by the particles /. The critical currents are 
determined by the intersection of the CVD and 2% of the flux flow voltage (dashed line). The optimal concentration of particles 
is at A^p = 4 000 corresponding to / = 0.23. (b) The dependences of the critical current Jc on the particle volume fraction / for 
a = 4^ and three magnetic fields. The optimal / increases with the magnetic field, (c) The dependences of Jc on / for different 
particle diameters at B = 0.1iifc2. The optimal particle diameter is in between 3^ and 4^. (d) Order parameter isosurfaces for 
a pinned vortex conhgurations with a = 4^, Np = 500, and B = 0.016//o2. (e) The field-induced vortex lines extracted from 
the same order parameter. The particles are shown as transparent spheres. Vortex lines outside particles are red, and inside 
particles are blue. 


as singularities in the order parameter. Even though 
the TDGL model does not provide a fully realistic de¬ 
scription of the dynamic properties of superconductors, 
it does describe accurately vortex-line flexibility, interac¬ 
tions between vortex lines, and interactions of vortices 
with pinning sites. It also allows for cuttings and recon¬ 
nections of vortex lines. Therefore, this model is per¬ 
fectly suited for the problem of the critical-current op¬ 
timization, for which a fully accurate description of the 
dynamics is not essential. The TDGL model has been 
proven to be very useful for exploring many properties 
of the vortex state [26-33]. However, only recently a 
meaningful exploration of the parameter space for suffi¬ 
ciently large three-dimensional (3D) superconductor has 
became possible [34] allowing us to address the problem 
of critical-current optimization. 

We use the TDGL model to explore vortex pinning 
by randomly-placed metallic spherical inclusions. Our 


objective is to find the optimal parameters for the pin¬ 
ning landscape to maximize the critical current. Clearly, 
when the particles occupy a small fraction of the total 
volume, the critical current grows as particle density in¬ 
creases. At some density, however, further increase of 
particle number will not improve current-carrying capac¬ 
ity [18] due to at least two factors: (i) The increasing 
mobility of the vortex lines due to jumping between the 
particles and (ii) the reduction the effective cross section 
for the supercurrent caused by inclusions. Therefore, it 
is important to find the size and density of particles that 
maximizes the critical current. This problem can not 
be accessed by simple approaches. We find optimal pa¬ 
rameters for different magnetic fields by systematically 
exploring the dependence of the critical current on size 
and concentration of particles. 

The dynamics of the order parameter 'ijj{r,t) is de- 
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scribed by the TDGL equation in the reduced form 

{dt + 

+ ^ + ( 1 ) 

j=x,y,z 

Here /r and A are the scalar and vector potentials. We 
used the in-plane coherence length ^ at the chosen tem¬ 
perature as the unit of length, meaning that rj^ = rjy = 1 
and rjz = I/7, were 7 is the anisotropy factor. We 
took 7 = 5 corresponding to YBCO. The function e(r) 
models pinning centers, e(r) = 1 in the bulk [35] and 
e(r) = —1 inside metallic inclusions. We used the approx¬ 
imation of large London penetration depth A in which 
the vector potential A is fixed by the external mag¬ 
netic field. Ay — Bx, and the magnetic field is mea¬ 
sured in units of the c-axis upper critical field at given 
temperature, Hc 2 = ‘f’o/(27r^^). The Langevin term 
C(r, t) describing thermal noise has the correlation func¬ 
tion (C*(r, t')) = T ( 5 (r—r') S(t—t'), where T is the 

reduced temperature in units of (Stt) and He is the 

thermodynamic field. The total electric current in units 
of Jo = c<i)o/(87r^^A^) (CGS) is given by 

J, = r,] {Imir (V, - - V,^x). (2) 

where the first term is the supercurrent and the second 
term gives the normal current. In these units the de¬ 
pairing current is Jdp = 2/(3-\/3 )Jq « 0.385Jo. We per¬ 
formed simulations with fixed current applied in the x di¬ 
rection and compute the average electric field E = — Vj,/! 
in the dynamic steady state. For the simulations, we 
developed a stable and efficient solver implemented for 
graphics processing units (GPUs) allowing for simula¬ 
tions of quite large 3D systems with sizes up to 500 co¬ 
herence lengths in all directions [34]. 

We systematically computed the current-voltage de¬ 
pendences (CVDs) for different particle sizes and den¬ 
sities. The simulated system size is 100^ x 100^ x 50^ 
with 256 X 256 x 128 mesh points. Figure 1(a) shows 
representative series of CVDs for different numbers of 
particles Ap with diameter of a = 4^ at magnetic field 
B — {)AHe 2 corresponding to 159 vortex lines in the sys¬ 
tem. These CVDs are obtained by stepwise decrease of 
the applied current with a simulation time of about 10® 
TDGL units between current steps. Typically, we did not 
observe significant history effects: CVDs computed with 
different starting currents and current steps are found to 
be very similar to each other. Only at small magnetic 
fields ^ 0.01Jlc2 CVDs become more noisy and slightly 
history-dependent. 

The pinning effectiveness of the particles is primarily 
determined by the volume fraction / occupied by them. 
For spatially separated particles the “nominal” volume 
fraction is /„ = irNpU^/{6LxLyLz). As in our case 
randomly-placed particles may overlap, the real volume 


fraction is somewhat smaller, f k, f^ — /^/2, where the 
correction term accounts for possible overlaps between 
pairs of neighboring spheres. For each number of parti¬ 
cles, the value of / is specified in parenthesis in Fig. 1(a). 
For low / the CVDs systematically shift to the right as 
particle density increases indicating an increase of the 
critical current. Above a certain density adding more 
particles starts to degrade the critical current. The opti¬ 
mal density corresponds to the volume fraction / = 0.23. 

The computed CVDs are used to evaluate the critical 
currents, Jc, which we define using the criterion A (Jo) = 
0.02pff Jo [dashed line in Fig. 1(a)]. The free flux-flow 
resistivity, ps, in our reduced units is ps = 1.689H. 
Figure 1(b) shows the dependences of the critical cur¬ 
rent on the nonsuperconducting volume fraction / for 
0 = 4^ and three magnetic fields. The optimal volume 
fraction slowly increases with the magnetic field, from 
~ 0.15 at H = 0.016iLc2 to ^ 0.23 at B = 0.1iLc2- At 
the lowest field, B = 0.016iLc2j the maximum current 
is Jc « 0.0383Jo corresponding to 10% of the depairing 
current. The movie in [36] demonstrates the evolution of 
vortex dynamics with increasing particle density for this 
magnetic field. 

Figure 1(c) presents the dependences of the critical cur¬ 
rents Jc on the nonsuperconducting volume fraction / for 
different particle diameters at B = 0.1iLc2- The highest 
critical current is realized at a = 3^ indicating the exis¬ 
tence of an optimal particle size. By applying parabolic 
fits to the numerical data for the geometry dependent 
Jc(o,/), we found the optimal particle size and volume 
fraction at several fields. Figure 2(a) shows the magnetic 
field dependences of the optimal particle size Oopt and vol- 



Figure 2. (a) The magnetic field dependence of the optimal 

diameter (left axis) and volume fraction (right axis), (b) The 
maximum critical current for optimal parameters for different 
magnetic fields in reduced units. The right axis shows this 
current normalized to the depairing current. 
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ume fraction /opt in the range 0.0167Jc2 < B < 0 . 2 Hc 2 - 
Within this range the optimal size monotonically de¬ 
creases with the increasing magnetic field from ^ 4.5^ to 
^ 2.5/. This indicates that the typical scale of disorder 
has to be comparable with the intervortex spacing which 
decreases as The optimal volume fraction has a 

weak, non-monotonic dependence on the field strength, 
but stays within the range of 17-22%. The magnetic-field 
dependence of the maximum critical current achieved for 
the optimal parameters is shown in Fig. 2(b). 

To characterize the structure of the pinned vortex 
states, we extracted the held-induced vortex lines from 
the order parameter [37] and performed a detailed anal¬ 
ysis of these conhgurations. Figure 1(d) shows a repre¬ 
sentative trapped vortex configuration for B = 0.016i7c2 
imaged by the order-parameter isosurfaces |'0(r)| = 0.1. 
Both particles and vortex lines can be seen as regions 
of suppressed order parameter. We found that the vor¬ 
tex arrangements typically are quite disordered, which is 
partly caused by the interaction of the vortices with ran¬ 
domly arranged particles and partly by incomplete equi¬ 
libration. Vortex lines traced from this order-parameter 
distribution are shown in Fig. 1(e). They are split into 
line segments located inside the particles and in super¬ 
conducting material, as illustrated by blue and red lines, 
respectively. 

We extracted several parameters characterizing trap¬ 
ped configurations: (i) The fraction of particles occu¬ 
pied by vortices, /fin, (h) the fraction of particles double- 
occupied by vortices, , (hi) the fraction of the total line 
length located outside particles, /free = ^outside/4otai, 

(iv) the average length of line segments trapped between 
neighboring particles, Lt, see inset in Fig. 3(c), and 

(v) the average particle-to-particle displacement in the 
direction of motion, ui = Uy, and in transverse direction 
Ut = Ux- These parameters are not independent. Indeed, 
the number of particles holding a given vortex line can be 
estimated as /fieeT^/Lt meaning that the total number 
of occupied pins is N^ftreeLz/Lt, where fVv is the total 
number of vortex lines. As fraction of pins /2 is holding 
two vortex lines, we can estimate the occupied fraction 
as /fill « Vv/free7^2/[fVp(l -f / 2 )Tt]- We checked that the 
extracted parameters satisfy this consistency condition. 

Particle fraction occupied by the vortex lines, /fin, nat¬ 
urally characterizes the efficiency of pin utilization. Fig¬ 
ure 3(a) shows the dependence of this parameter on the 
nonsuperconducting volume fraction / for a = 4/ and 
three magnetic helds. As expected, this parameter in¬ 
creases with the magnetic field and decreases with the 
number of particles. For B = 0.1i7c2 almost all particles 
are occupied. We also found that for this field typically 
3-5% of particles hold two vortex lines without any sys¬ 
tematic dependence on the particle density. For smaller 
fields the number of double-occupied particles is negli¬ 
gible. An noteworthy feature is that for small fields a 
significant fraction of particles remains unoccupied even 



Figure 3. Evolution of parameters characterizing pinned 
vortex-line configurations with increasing nonsuperconduct¬ 
ing volume fraction / for a = 4/ and three magnetic fields. 
The plot (a) shows the fraction of particles occupied by vortex 
lines, /fill. The plot (b) presents the length fraction of vor¬ 
tex segments outside the particles, /free- Dashed line shows 
the volume fraction occupied by superconducting material, 
1 —/. The plot (c) shows the average geometrical parameters 
characterizing pinned line segments as illustrated in the in¬ 
set, the segment length Lt (left axis) and typical pin-to-pin 
line displacements in the direction of motion, ui, and in the 
transverse direction, Ut (right axis). The parameters Lt and 
Ml are defined in the inset. In all plots arrows mark locations 
of the maximum critical current. 


for very low particle densities. 

The efficiency of the vortex trapping by the particles 
can be characterized by the free-segment length fraction 
of the vortex lines, /free, plotted in Fig. 3(b). A natural 
upper limit for /free is the volume fraction occupied by 
the superconducting material, 1 — /, shown by the dashed 
line. This limit would be realized if there were no corre¬ 
lations between particles and vortices. We can see that 
for pinned configurations /free is significantly below this 
limit since particles trap vortices. Notably, /free drops 
below 50% when particles occupy only 5% of the volume 
for B = 0.016i7c2- The free-segment fraction monotoni¬ 
cally increases with the magnetic field, because the vor¬ 
tex lines compete for the particles. We observe that the 
maximum critical current is realized for /free = 21-25% 
which only weakly depends on the magnetic field. 
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For strong pinning sites the vortex lines split into 
finite-size segments hanging in between neighboring 
sites [13, 14, 23]. Figure 3(c) shows the behavior of the 
average length parameters in units of ^ characterizing 
these free-line segments, the segment length Lt and pin¬ 
to-pin line displacements in the direction of motion, ui, 
and in the transverse direction, ut- The length Lt rapidly 
increases with decreasing particle density and with in¬ 
creasing magnetic field. While the displacements t also 
increase with decreasing /, they only weakly depend on 
the magnetic field. As expected, for small particle den¬ 
sities the vortices stretch between pinning sites preferen¬ 
tially in the direction of motion [13, 14, 23] meaning that 
ui is larger than ut- These parameters, however, become 
almost identical for / > 0.15. At the particle density 
corresponding to the maximum of the critical current all 
three parameters, Lt, u\, and itt, are close to ^ for both 
magnetic fields. 

In summary, we investigated vortex pinning by ran¬ 
domly distributed metallic inclusions and found optimal 
parameters for highest critical current. We also analyzed 
statistical properties of pinned vortex arrays and revealed 
several nontrivial structural properties of the optimally- 
pinned states. Our general observation is that there is no 
universal optimal pinning configuration for all magnetic 
fields. Thus, for best performance in a given applica¬ 
tion, pinning landscapes should be designed taking into 
account relevant magnetic fields. 
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